Interaction of cosmic background neutrinos with matter of periodic structure 
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We study coherent interaction of cosmic background neutrinos(CBNs) with matter of periodic 
structure. The mixing and small masses of neutrinos discovered in neutrino oscillation experiments 
indicate that CBNs which have very low energy today should be in mass states and can transform 
from one mass state to another in interaction with electrons in matter. We show that in a coherent 
scattering process a periodic matter structure designed to match the scale of the mass square 
difference of neutrinos can enhance the conversion of CBNs from one mass state to another. Energy 
of CBNs can be released in this scattering process and momentum transfer from CBNs to electrons 
in target matter can be obtained. 

PACS numbers: 13.15.+g 



CBN is one of the major predictions of the big- 
bang cosmology. The big-bang cosmology predicts 
that the temperature of CBNs, if they are relativis- 
tic today, is around 1.96 K |l| (~ 10 -4 eV) in the 
present universe. The average number density of 
CBNs per species in the present time is predicted 
to be 



56 cm 
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The detection of CBNs is extremely difficult due 
to the extremely low energy and the small density 
of CBNs in the present universe. 

In this article we study the scattering of mas- 
sive CBNs with matter of periodic structure. We 
point out that a periodic structure of matter, when 
matching with the scale of the mass square differ- 
ence of massive CBNs, can help to convert CBN of 
one mass state to another and enhance the prob- 
ability of CBNs scattering with matter. Net mo- 
mentum transfer from CBNs to target matter can 
be achieved in this coherent scattering process. 

In the following we first have a brief review of 
present knowledge of neutrino masses and mixing. 
Then we study the scattering and conversion of 
CBNs in a target matter of periodic profile. We 
discuss momentum transfer from the wind of CBNs 
to target matter. Finally we conclude. 

Neutrino oscillation experiments have shown 
that neutrinos have very small masses and flavor 
mixing. Two mass square differences of neutrinos 
have been measured in oscillation experiments [2| : 
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Flavor mixing of neutrinos is expressed as 



where v\ is neutrino state in flavor base and Vi 
neutrino state in mass base. Neutrino mixing ma- 
trix U can be parameterized using three mixing 
angles #12,23, 13 and a CP violating phase 0. In 
solar, atmospheric and long baseline neutrino os- 
cillation experiments two mixing angles have been 
measured: sin 2 26 u « 0.86, sin 2 6> 23 ps 0.50 
. Recent observation of oscillation of reactor anti- 
neutrinos in Daya Bay experiment [3j, confirmed 
by RENO experiment [J] , shows that sin 2 2(9i3 ps 
0.092. According to all these measurements of neu- 
trino oscillation we can conclude that no element 
of U is zero. 

Neutrino masses are also measured in /3 decay 
experiment and in cosmological observations, ft 
decay experiment gives an upper bound on the 
mass of electron anti-neutrino [2J 



< 2.0 eV 
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The observations of the anisotropy of Cosmic Mi- 
crowave Background Radiation and the Sloan Dig- 
ital Sky Survey give a constraint on the total mass 
of neutrinos [2| 



E 



< 0.8 eV 
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Taken all these measurements of neutrino masses 
into account, possible patterns of neutrino masses 
are: Normal Hierarchy (NH) with m\ <C 1112 < W3, 
v/Am^ fa 0.9 x 10~ 2 eV and m 3 ps 
0.05 eV; Inverted Hierarchy with 
« 0.05 eV; 
m 2 ss 7713 < 0.3 
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m 3 <C mi < ni2, mx t 2 ~ y/\Am 
Quasi-Degeneracy(QD) with m,\ 
eV, mi > VI A ' 
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Massive CBNs with m > 10" 4 eV should be 
non-relativistic in the present time and should be 
in mass states. Massive CBNs with low veloc- 
ity should also be clustered in galactic halos or in 
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cluster halos. For different neutrino mass patterns 
massive neutrinos arc different and the contents of 
non-relativistic CBNs are different too. According 
to the above discussion of neutrino mass patterns 
we can figure out that the CBNs which are non- 
relativistic today include 1/2,3 for NH, v\^ for IH 
and all 1/1,2,3 for QD. In the following we will con- 
sider interaction of massive CBNs with matter. 

The interaction of neutrinos relevant to our anal- 
ysis is the interaction which is non-universal in 
neutrino flavors. Neglecting radiative corrections 
neutral current interaction of neutrinos with mat- 
ter is universal in flavors and is irrelevant to later 
analysis. The relevant interaction is given by the 
charged current interaction of neutrino with elec- 
tron: 




FIG. 1: Scattering of neutrino with matter of periodic 
structure. 
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where and v e are the fields of electron and elec- 
tron neutrino with left-chirality. Gf is the Fermi 
constant. Using ([7]) one can find that in an un- 
polarized target of matter at rest the coherent in- 
teraction of neutrino with matter gives a potential 
term to electron neutrino: 



AC = -V2G F N e v e ~f°v e 



(8) 



where N e is the number density of electron in mat- 
ter. (JSJ describes coherent scattering of neutrino 
with matter in which neutrino coherently scatters 
with many electrons in matter. Momentum trans- 
fer from neutrino, if not zero, is distributed to very 
large numbers of electrons participating actively in 
the scattering process and momentum transfer to 
a single electron can be taken zero. This case is 
exactly what we study in later discussion. 
In the mass base (TSl) can be rewritten as 



AC 



-V e U* e] U el Vrfvi, 
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where i,j = 1,2,3 and V e = V2GfN c . According 
to ([5]) , neutrino in one mass state z/j can transform 
to another mass state Vj in interaction with elec- 
trons in matter. Consider such a transition — > Vj 
in a target of matter. For a uniform incident flux of 
neutrino Vi the cross section for the Vi — Vj (pi — Vj) 
conversion in matter is 
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where Ei and Ej are the energies of the initial 
Vi and final Vj respectively, Vi the velocity of the 
initial i>{ relative to the target, fcj and kj the 
initial and final momenta, is the volume of 
the target. \M\ 2 is the matrix element squared: 



\M\ 2 = k®k® + ki ■ kj for unpolarized Dirac type 
neutrinos. In later discussion we will concentrate 
on massive Vi CBNs which are non-relativistic and 
are clustered in present time. Since the direction 
of motion is changed in clustering process but the 
spin of neutrino is not, the clustered CBNs can be 
considered mixed with left and right helicities. So 
we can take <C kf and use \M\ 2 — fc°/c° in 
this case. For Majorana neutrino the neutrino and 
anti-neutrino are identical and \M\ 2 is replaced by 
\M\ 2 = 2(kikj + ki ■ kj — niimj). It is velocity sup- 
pressed if Vi and Vj are both non-relativistic. For 
relativistic Vj and non-relativistic Vi one can use 
the approximation \M\ 2 = 2EiEj. For simplicity 
we will concentrate on Dirac type neutrino and use 
\M\ 2 = EiEj in later discussion. 

We consider a target of matter which is constant 
in x and y directions and periodic in z direction, as 
shown in Fig. [TJ The potential term in such kind 
of matter profile satisfies 



V e (z + d) = V e (z), 



dx 



dVe_ 

dy 



= 0, 
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where d is the period of the matter profile. A gen- 
eral potential term of such kind periodic structure 
can be expressed using Fourier transformation as 



V e (x) = J2 V n e l 



(12) 



where n is an integer, q n = q n z and q n — 2nir/d. 
V n satisfies: V* = V- n . 

If the path length of neutrino in target is con- 
stant the cross-section can be expressed as a = Sp 
where p is the probability of Vi — Vj conversion 
and 5* is the geometric cross section of the target 
of matter. Implementing (TTSl into (TIU)) and in- 
tegrating over space coordinates we find that the 
conversion probability is 



\k z A\M\ 2 

^=-4^k lVnLzU:jUe 



Asm 2 (A n L z 
(A n L z ) 2 
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where 

A„ = k\ - k] - q n , (14) 

and L z is the length of target in z direction. For 
T>i — Vj conversion the probability is the same. 
We can see that the probability is proportional to 
|U„/A„| 2 if \A n L z \ > 1. If \A n L z \ « 1 the con- 
version is resonantly enhanced and the probability 
is proportional to |VT„ J L Z | 2 . At the resonant point 
of V{ — v,j conversion we can find 

kf = fcj , k^ = kj, fcj kj — q n = (15) 

This condition of resonant conversion is expressed 
in short as 

h - kj - q n = 0. (16) 

When the resonant conversion happens the cross- 
section and the probability p n are proportional to 
\V n | 2 of one particular n and other V n i-t n effectively 
contribute zero. 

Using energy conservation we can find that 

&| = + $ ( 17 ) 

Using kf — kj, k\ = fcj and (|17l) we can find that 

= V^f - m] + {klf (18) 

We will concentrate on massive CBNs which are 
non-relativistic in the present universe. The veloc- 
ity of massive CBNs at the position of solar system 
depends on the clustering properties of CBNs. If 
CBNs are in virial equilibrium in the local galac- 
tic halo the velocity of these CBNs is ~ 200 km/s 
~ 10 _3 c where c is the speed of light. If CBNs 
are in virial equilibrium in the local cluster halo 
the velocity of CBNs is a bit larger and can reach 
~ 1000 km/s ~ 10~ 2 c. According to © we can 
get rrii < 0.2 — 0.3 eV. So we can conclude that 
= [rriiVil < 10 -3 eV and kf < |Am| 132 |. So 

using (USD we get |/cf| « ^Am?-. \k?\/Ej in (T3) 

can be approximated as sa yjl — rrVj/m?. 

Apparently the i>i — Vj conversion can not happen 
if Am 2 , < 0. Form 2 — m 2 > 0, Vi — Vj conversion is 
allowed and it means that non-relativistic z/j con- 
verts to Vj and releases part of its rest energy to 
kinetic energy of Vj . 

The transferred momentum from neutrino to the 
target is fcj — kj. Using (fTo) we can find that for 
resonant i>i — Vj conversion the momentum transfer 
is q n . Using (ITSl) we can find that \q n \ — \kj—kf\ « 

y^Am 2 " and q n m ±^J AmfjZ. q n can be positive 

or negative which correspond to cases that Vj is 
reflected or is refracted by the target matter. The 
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FIG. 2: Asymmetry of p'±x versus for vi — V\ tran- 
sition of NH: mi = ^/Am^. k s L z =5. k s — 77i2|# s |- 

net momentum transfer per unit time from CBNs 
is 

P= \q±\\S m J dvi (p+i— p-i)vif(Vi), (19) 

where / is the velocity distribution of local 
CBNs, rii the number density of local Vi and Pj 
CBNs, S the geometric cross section of target mat- 
ter. P can be positive or negative which corre- 
spond to cases that the momentum transfer to 
target is of positive or negative z direction. Ap- 
parently there will be no net momentum transfer 
from CBNs to target detector if the probabilities of 
CBNs being refracted or being reflected by target 
matter are equal. A periodic structure of target 
matter can make these two probabilities differ sig- 
nificantly and net momentum transfer from CBNs 
to target can be obtained. 

To illustrate that net momentum transfer can be 
achieved we consider, as an example, the case that 
q+i t -i dominate the i>i — Vj conversion. In such 
a case the period of the target detector should be 

arranged to satisfy ^ ~ Amfj . When neutrino 

is refracted by the detector fc| = Am 2 ^ + (fcf) 2 
and n = — 1. We get 

A_j = fc? + y " V A <" + ^) 2 - ( 2 °) 
When neutrino is reflected by the detector fcj = 
— ^/Am?- + (fcf) 2 and n = +1. We get 

A+i = *f - y + \J Am iJ + ( fc *) 2 - ( 21 ) 

We can see that p+i — p~i would not be zero and 
net momentum transfer would be obtained if the 
detector is arranged in such a way that |A + i_L^| < 
1 < |A_iL,| or |A_iL,| < 1 < \A +1 L Z \. 

To achieve net momentum transfer from neu- 
trino background we note that the solar system is 
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FIG. 3: Ap' x k s L z v ersus v z a for ^2 — ^1 transition and 
QD: m,2 = y^AmfTl- Three lines are for k s L z — 80, 40, 
10, 5 and 2 separately. k 3 — m,2\v 3 \ 



FIG. 4: Asymmetry of p±i versus |g±i| for ^3 — i^i 
transition for QD: m-j = 3 x Am^ [. fc 3 L z = 10. 
fcf = fc s = m3|v s |. Values of q±i| for resonant conver- 
sion varies for different Amli. 



moving in the local galactic halo or cluster halo. 
The momentum of CBNs in the rest frame of solar 
system can be written as 



k s ~t~ Aki : 



(22) 



where k s is the momentum caused by motion of 
solar system relative to CBN halo. If the structure 
of detector is arranged such that \q±i\ — = 

rn^A + S we can find that 



1 (k z ) 2 

A_ x = 2k z s + Afcf + 6 - - \ * ' , (23) 



A +1 = Akl-5 + l- {kf)2 



(24) 



where 5 is a possible small mismatch between =?■ 
and k z s + ^Am 2 l3 . Since \k i \/JArr^ j < lO^-lO" 2 

as observed in previous discussions the last terms 
in (H3J) and (J2H) can be neglected. If \2k z L z \ > 1 
it's easy to see that a difference between p +1 and 
p_i can be achieved. Apparently if k z s = two 
probabilities should be equal. In Fig. [2] we can see 
the asymmetry of p±\ clearly when v z approaches 
250 km/s. 

In FigJ3]we compute Ap 1 : 



where \M | 2 rj EiEj has been used. As an estima- 
tion we use the Maxwellian distribution 



/(«) = 



1 



= -(^-<) 2 /a 2 



(27) 



In galactic halo v z can reach |u s | = 250 km/s when 
2 direction follows the direction of motion of the 
solar system in Milky Way. a is the velocity dis- 
persion of CBNs which is taken as a = 150 km/s 
in galactic halo. In computing Fig. [J] we have 
used 2it/d = m2\v s \ + y/ Am^ and 8 = 0. We can 
see in Fig. [5] that the asymmetry of p' ±1 disap- 
pears when the v z approaches zero and it's maxi- 
mal when v z = \v s \. 

Ap' does not reach 1 partly due to the cancella- 
tion of p'±x- For |A: S L Z | < 1 the cancellation can be 
significant. For |fc s L z | ~ 1 we can have Ap 1 of or- 
der one. The other factor which reduces Ap 1 is that 
for large L z only parts of CBNs with \Ak z L z \ < 1 
contribute to the resonant conversion. The larger 
L z is, the smaller the fraction of CBNs contribut- 
ing to resonant conversion. \Ap'k s L z \ actually in- 
creases as k s L z increases and it approaches a fixed 
value for k s L z 3> 1. In Fig. [3] we can see this 
effect. We can find that for k s L z = 80 and 40 
the two curves of Ap'k s L z converge. In this case 
Ap = p + i — p_i approaches to a value which is 
approximately proportional to 3L z /k s 

We note that to achieve maximal Ap' sufficiently 

good matching between \q±i\ and ^/ Am^ ±irii\v s 



d"i /«) (p'+i-p'-i), (25) is needed. Unfortunately we do not have very pre 

• _ 1 _r A„„2 /FTl m.j_ i j j.. 



where p' n = s'm 2 (A n L z )/(A n L z ) 2 . The net mo- 
mentum transfer per unit time from CBNs, P, can 
be expressed using Ap' as 



TYl 

P = \q±i\Sn iC \ll- -^\V ±1 L z U* ej U ei \ 2 Ap', (26) 

TYl: 



cise knowledge of Amy. With ([3]) it's hard to 
achieve a matching with precision to one of a hun- 
dred or one of a thousand. To overcome this prob- 
lem one can use a number of detectors which makes 
a scan of the range of Amfj . For example one can 
use a hundred copies of detectors with identical L z 
and slightly different |<7±i|. The values of |g±i| are 
evenly distributed in the uncertain range of Am| . 
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FIG. 5: Asymmetry of p'±x versus |g±i| for — v\ 
transition. k a L z = 10. k\ = k 3 — m2\v s \- NH: m,2 = 

3 x 



s jAm\ x ; QD1: m 2 = 2 x y'lAmfJ; QD2: 
V\ Am 3i\- 
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In Fig. [4] we give an example for 1/3 — v\ conver- 
sion and QD: TO3 = 3. * ^/|Am|jJ- v "i i s considered 
heavier than v\ in this example. The range of |<7±i| 
in this figure corresponds to lAm^l in the uncer- 
tain range [2.08, 2.68] x 10" 3 eV 2 . We see that the 
values of |<7±i| which give resonant conversion are 
different for different values of | Am^ | and for each 
lAmfj-J there are 6 — 8 values of |g±i| of total one 
hundred which give resonant enhancement. Posi- 
tive or negative Ap' correspond to the cases that 
the momentum transfer to target is of the positive 
or negative z direction. We note that by carefully 
adjusting the period of target matter the momen- 
tum transfer from CBNs can be of the positive or 
negative direction of the CBN wind. 

The number of |<?±i| which gives resonant con- 
version depends on the neutrino mass pattern. For 
QD mass pattern the neutrino mass is larger than 
that of NH and IH. The initial momentum of CBNs 
are larger too. So it's easier for QD to achieve a 

matching of \q±i\ and Amf^ ± k s to the preci- 
sion of k s . In Fig. Uwe can see this clearly in 
i>2 — v\ conversion. The range of |<?±i| in this fig- 
ure corresponds to Ato 2 ^ in the uncertain range 
[6.90,8.17] x 10- 5 eV 2 . We see that for QD the 
resonant region is broad. But for NH the resonant 
conversion happens in a narrow region of |?±i|- 

For matter on Earth we find that V e ~ 10~ 13 
eV. A detector can be designed to have V±i of the 
same order of magnitude of V e . For L z = 1 cm we 
have |V±iL 2 | ~ 10~ 8 . When sufficiently large Ap' 
is achieved we can have a rough estimate of the net 
momentum transfer per unit time from CBNs: 




(28) 



This momentum transfer to target matter gives 
rise to a mechanical force exerted on the target 
detector. It's easy to estimate that the acceler- 
ation due to momentum transfer (|2"8")l can reach 
10 -27 — 10 -28 cm s~ 2 . As can be seen the net 
momentum transfer to matter is very small. The 
mechanical force exerted on the target is also very 
small. It's far less than that can be detected in 
modern technology. Moreover, it might be much 
smaller than the possible mechanical force exerted 
by dark matter [5|. Detecting this momentum 
transfer from CBNs using mechanical force is not 
possible at the moment. We should think about 
other mechanisms to detect the momentum trans- 
fer from the wind of CBNs. We note that the mech- 
anism considered in this article is only sensitive to 
the charged interaction of neutrino with electrons 
which is non-universal in neutrino flavors. The mo- 
mentum is transferred from CBNs to electrons in 
matter through coherent scattering process. Con- 
sidering electrons in conductor or superconductor 
might lead to a better way to detect momentum 
transfer from the wind of CBNs. The research of 
this topic is out of the scope of the present article. 

We note that the mechanism discussed in this 
article make uses of the fact that massive neutri- 
nos can convert from one mass state to another in 
interaction with electrons in matter. Momentum 
is transferred from CBNs to electrons in coherent 
scattering process with target matter. This is dif- 
ferent from previous works 043 which considered 
the coherent scattering of neutrino with nuclei in 
matter. Another difference from previous works is 
that the momentum transfer from CBNs to target 
matter discussed in this article can be of positive or 
negative direction of the CBN wind, depending on 
the period of the target matter. Momentum recoil 
given by CBN wind to target matter, discussed in 
p-uj , is always of the same direction of the CBN 
wind. 

We note that if n% in galaxy is much larger than 
the average density fl} the event rate can be en- 
hanced. However, due to Pauli blocking Q it's 
difficult for the density of CBNs in our galaxy to 
be much larger than the value in ([1]) unless the 
neutrino mass reach ~ 0.5 eV @. Numerical sim- 
ulation of clustering of CBNs does not support m 
in galaxy much larger than the average value ei- 
ther 10]. If considering CBNs clustered in local 
cluster halo the number density is allowed to be 
much larger and the signal of CBNs is larger. The 
mechanism considered in this article can also be 
applied to eV scale sterile neutrino or keV scale 
sterile neutrino dark matter when the period of 
the target matter is designed to be the scale of //m 
or nm. 

In conclusion we have studied the coherent scat- 
tering of CBNs with a detector of periodic mat- 
ter structure. Massive CBNs which are non- 
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relativistic today can convert from one mass state 
to another mass state in interaction with electrons 
in matter. Energy of neutrino is released in this 
scattering process and momentum can be trans- 
ferred from CBNs to target matter. We show that 
a periodic structure of matter can enhance the 
scattering probability when the period is matched 
to the scale of the mass square difference of neu- 
trinos. A good arrangement of the periodic struc- 
ture can also select the CBNs to be reflected or 



be refracted by the target matter and lead to net 
momentum transfer to the target matter from the 
wind of CBNs. If a smart way to detect this small 
momentum transfer can be found the result found 
in this article might be useful for designing a real- 
istic detector for detecting CBNs in laboratory. 
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